Chaos in Superstring Cosmology 



Thibault Damour 

Institut des Hautes Etudes Scientifiques, 35, route de Chartres, F-91440 Bures-sur-Yvette, France 

and 

Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030, USA 

Marc Henncaux 

Physique Theorique et Mathematique, Universite Libre de Bruxelles, CP. 231, B-1050, Bruxelles, Belgium 

and 

Centro de Estudios Cientificos, Casilla 1469, Valdivia, Chile 

It is shown that the general solution near a spacelike singularity of the Einstein-dilaton-p-form field 
equations relevant to superstring theories and M-theory exhibits an oscillatory behaviour of the 
Belinskii-Khalatnikov-Lifshitz type. String dualities play a significant role in the analysis. 



An outstanding result in theoretical cosmology has 
been the discovery by Bclinskii, Khalatnikov and Lifshitz 
(BKL) that the generic solution of the four-dimensional 
Einstein's vacuum equations near a cosmological singu- 
larity exhibits a never ending oscillatory behaviour |Q. 
The oscillatory approach toward the singularity has the 
character of a random process, whose chaotic nature has 
been intensively studied However, two results cast 
a doubt on the physical applicability, to our universe, of 
the BKL picture. First, it was surprisingly found that 
the chaotic BKL oscillatory behaviour disappears from 
the generic solution of the vacuum Einstein equations in 
spacetime dimension D > 11 and is replaced by a mono- 
tonic Kasner-like power-law behaviour || . Second, it was 
proved that the generic solution of the four-dimensional 
Einstein-scalar equations also exhibits a non-oscillatory, 
power-law behaviour Q , || . 

Superstring theory || suggests that the massless 
(bosonic) degrees of freedom which can be generically 
excited near a cosmological singularity correspond to a 
high-dimension (D = 10 or 11) Kaluza-Klein-type model 
containing, in addition to Einstein's £>-dimensional grav- 
ity, several other fields (scalars, vectors and/or forms). 
In view of the results quoted above, it is a priori un- 
clear whether the full field content of superstring the- 
ory will imply, as generic cosmological solution, a chaotic 
BKL-like behaviour, or a monotonic Kasner-like one. In 
this letter we report the result that the massless bosonic 
content of all superstring models (D = 10 IIA, IIB, I, 
hetE, hetso), as well as of M-theory (D = 11 supergrav- 
ity), generically implies a chaotic BKL-like oscillatory 
behaviour near a cosmological singularity. [Our analy- 
sis applies at scales large enough to excite all Kaluza- 
Klein-type modes, but small enough to be able to ne- 



glect the stringy and non-perturbative massive states.] 
It is the presence of various form fields (e.g. the three 
form in SUGRAn) which provides the crucial source of 
this generic oscillatory behaviour. 

Let us consider a model of the general form 

S = J d D x y /g[R(g)-d fl ipd^ 

Here, the spacetime dimension D is left unspecified. We 
work (as a convenient common formulation) in the Ein- 
stein conformal frame, and we normalize the kinetic term 
of the "dilaton" ip with a weight 1 with respect to the 
Ricci scalar. The integer p > labels the various p-forms 
A p = A^...^ present in the theory, with field strengths 
F p+1 = dA p , i.e. w ... Mp = d tlQ Ap 1 ...p p ±p permu- 
tations. The real parameter A p plays the crucial role of 
measuring the strength of the coupling of the dilaton to 
the p-form A p (in the Einstein frame). When p = 0, we 
assume that Ao ^ (this is the case in type IIB where 
there is a second scalar). The Einstein metric g^ is used 
to lower or raise all indices in Eq. (Q) (g = — det g^ u ). 
The model (|l|) is, as it reads, not quite general enough 
to represent in detail all the superstring actions. Indeed, 
it lacks additional terms involving possible couplings be- 
tween the form fields (e.g. Yang-Mills couplings for p = 1 
multiplets, Chern-Simons terms, {dCi — Co d £?2) 2 -type 
terms in type IIB). However, we have verified in all rel- 
evant cases that these additional terms do not qualita- 
tively modify the BKL behaviour to be discussed below. 
On the other hand, in the case of M-theory (i.e. D = 11 
SUGRA) the dilaton ip is absent, and one must cancell 
its contributions to the dynamics. 

The leading Kasner-like approximation to the solution 
of the field equations for g^ v and ip derived from ([!]) is, 
as usual W, 
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g pv dx" dx v ~ ~dt 2 + t2n(x) i^f > 



p ~ p v (x) In t + -0 (a;) . 
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where d = D — 1 denotes the spatial dimension and where 
lo 1 (x) = (x) dxi is a time-independent d-bein. The 
spatially dependent Kasner exponents Pi (x), p v (x) must 
satisfy the famous Kasner constraints (modified by the 
presence of the dilaton) : 



Q{p) = p 
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The set of parameters satisfying Eqs. (|J) is topologically 
a (d-l)-dimensional sphere: the "Kasner sphere" . When 



the dilaton is absent, one must set p v to zero in Eq.(3a). 
In that case the dimension of the Kasner sphere is d— 2 = 
D — 3. 

The approximate solution Eqs. (^) is obtained by ne- 
glecting in the field equations for <? M „ and tp: (i) the effect 
of the spatial derivatives of <?p„ and ip, and (ii) the con- 
tributions of the various p-form fields A p . The condition 
for the "stability" of the solution (§), i.e. for the absence 
of BKL oscillations at t — > 0, is that the inclusion in 
the field equations of the discarded contributions (i) and 
(ii) (computed within the assumption (||)) be fraction- 
ally negligible as t — ► 0. As usual, the fractional effect 
of the spatial derivatives of (p is found to be negligible, 
while the fractional effect (with respect to the leading 
terms, which are cx t~ 2 ) of the spatial derivatives of the 
metric, i.e. the quantities t 2 JJ*. (where iZ* denotes the 
e?-dimensional Ricci tensor) contains, as only "dangerous 
terms" when t — ► a sum of terms cx t 29ijk , where the 
gravitational exponents g^k (i ^ j, i 7^ k, j =/= k) are the 
following combinations of the Kasner exponents || 



9ijk (p) =2pi + 



^2 Pi = 1 + Pi - Pj - Pk ■ (4) 



The "gravitational" stability condition is that all the ex- 
ponents gijk (p) be positive. In the presence of form fields 
A p there are additional stability conditions related to the 
contributions of the form fields to the Einstein-dilaton 
equations. They are obtained by solving, a la BKL, the p- 
form field equations in the background (||) and then esti- 
mating the corresponding "dangerous" terms in the Ein- 
stein field equations. When neglecting the spatial deriva- 
tives in the Maxwell equations in first-order form d F = 
and 5 (e Ap v F) = 0, where S = * d * is the (Hodge) dual of 
the Cartan differential d and F p+ \ — dA p , one finds that 
both the "electric" components v / ge Ap v F 0n - lp , and the 



"magnetic" components Fj 1 ...j +1 , are constant in time. 
Combining this information with the approximate results 
(||) one can estimate the fractional effect of the p-form 
contributions in the right-hand-side of the g^- and ip- 
field equations, i.e. the quantities t 2 T? A -. Q and t 2 
where T^. v denotes the stress-energy tensor of the p- 
form. [As usual Q] the mixed terms T?^, which enter 
the momentum constraints play a rather different role 
and do not need to be explicitly considered.] Finally, 
one gets as "dangerous" terms when t — > a sum of 

2e (p) 

"electric" contributions cx t 11 and of "magnetic" 

ones cx t ^••■^-p-i . Here, the electric exponents e? 4 
(where all the indices i n are different) are defined as 



C (P) 



Pi!+Pl2^ 1" Pi 



while the magnetic exponents j d i (where all the 
indices j n are different) are 



2 ^pPcp 



(5) 



U- p -i (f) = Pn +Ph+'"+ Pm-p-i + 2 X vPf ■ ( 6 ) 



To each p-form is associated a (duality-invariant) double 
family of "stability" exponents e^ p \ The "electric" 
(respectively "magnetic") stability condition is that all 
the exponents e^ (respectively, b^) be positive. This 
result generalizes the results of on the effect of vector 
fields in D = 4. 

The main result reported here is that, for all super- 
string models, there exists no open region of the Kas- 
ner sphere where all the stability exponents g{p), e(p), 
b(p) are strictly positive. To define the set of stability 
conditions for the various superstring models, let us re- 
view their field content and give the values of the crucial 
dilaton couplings X p . The simplest case is the massless 
bosonic sector of M-theory, i.e of SUGRA in D = 11. In 
that case, there is a 3-form and no dilaton. The parame- 
ters p* f , a = 1, . . . , 10, run over the 8-dimensional sphere 
S S M defined by £ a (p™) 2 = 1 = £ a V% ■ Th e presence 
of a 3-form A3 uncoupled to any dilaton implies that the 
electric and magnetic stability exponents are respectively 
given by (||) and (||) with p = 3, X p — and d — 10, i.e., 

e ai a 2 a 3 — P ai +Pa 2 + Pa 3 ana °ai...a e — Pat + + Pa 6 ■ 

The D — 10 type IIA string theory involves, besides 
g pv and a dilaton p = $/\/2 (with g s — e* being 
the string coupling) a 1-form, a 2-form and a 3-form. 
The (Einstein-frame) dilaton coupling parameters of the 
forms are \± = 3^2/2, \£ = -V2 and Ag 4 = y/2/2, 
respectively. Besides the dilaton Kasner exponent p^, 
there are nine metric exponents pf, i = 1, . . . , 9. They 
run over defined by Eqs. (||). 

The D — 10 type IIB string theory involves (besides 
9nu)- two scalars: the dilaton p = ^j\/2 and the R — R 
0-form C , two 2-forms B 2 {NS-NS) and C 2 (R-R), and 
one "self-dual" R— R 4- form C4. The dilaton coupling 
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strengths of the forms are: A^ Q 



2^2, A| 2 - -V2, 
A^ 2 = +V2 and A^? 4 = 0. [A^ refers to the more compli- 
cated mixed coupling e* {dC 2 — Cq dB 2 ) 2 }- The Kasner 
exponents p^ , (i = 1, . . . , 9) run over S% defined by 
Eqs.®. 

The D = 10 type I string theory involves (besides g^ 
and if): an SO (32) vector potential, and a 2-form. The 
dilaton couplings are A{ = \/2/2 and X 2 = + \/2f]. The 
Kasner exponents p 1 ^, p\ (i — 1, ... ,9) run, as for IIA 
and IIB, on the Sf defined by Eqs. (§|). 

Finally, the D — 10 heterotic string theories involve 
(besides g^ and ip): an SO (32) or E% x E% vector po- 
tential, and a 2-form. Their respective (Einstein-frame) 
dilaton couplings are: \\ = — y/2/2, Af = —s/2. The 
Kasner sphere for p^, p^ (i = 1, . . . , 9) is the same as 
for IIA, IIB or I. 

Let us denote for each given theory "th" (where 
th = At, A, B, I, h labels the theory) the full (finite) se- 
quence of stability exponents as u»j h (p), where J labels 
all the possible exponents within each theory. E.g., when 
th = At the label J takes 690 values corresponding to 
the set {wf} = {s^ 7 , eaxaias, % e }- The condition 
of "Kasner stability" of each theory is that there exist 
an open region of the corresponding Kasner sphere S® h 
where w t ] 1 (p) > for all the labels J. However, we have 
proven that, for all theories, mijw t j 1 (p) is strictly nega- 
tive for all values of p € Sf h , except at a finite number of 
isolated points where it vanishes. 

Let us first consider M-theory. We have proven a 
stronger result, namely, that the electric stability condi- 
tions alone are never fulfilled. If, at any point on S^, we 



order the Kasner exponents as pf < p^ < 
the most stringent electric stability criterion involves 
/o (p) = pf 1 +p^ +P3 1 ■ To show that this function is non- 
positive on the cell pf 1 < ■ ■ • < pfg of the Kasner sphere, 
we maximize it subject to the constraints (Pa 1 ) 2 = 1j 
S Pa = 1- These constraints can be taken into ac- 
count by introducing two Lagrange multipliers. After a 
straightforward (but rather long) exhaustive analysis, we 
have found that / ( 5 nax — 0, this maximum being reached 
only at pi = ■ ■ ■ = p g = 0, p w = 1. 

To deal with the type IIA theory, we use the fact that 
IIA is the Kaluza-Klein (KK) reduction of M on a circle. 
This fact dictates the link between the field variables of 
the two models. If we label by the letter y the compacti- 
fied dimension this link implies the following relation be- 
tween the (Einstein-frame) Kasner exponents of the two 
theories (i = 1, . . . , 9) 
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8 P f +Pl 



(7) 



*Note the misprint in Eq. (12.1.34b) of 
author's web page. 



corrected on the 



Forgetting about this Kaluza-Klein motivation^ we can 
consider that Eqs. (fy define a one-to-one map %am from 
Sff to S%: p£ — tam {pff)- Using this map, we have 
then shown that the complete set of IIA stability con- 
ditions is logically equivalent to the complete set of M 
stability conditions. The instability of the Kasner be- 
haviour of M-theory proven above then implies that the 
Kasner behaviour of IIA is also unstable. 

To deal with the type IIB theory, we use the fact 
that IIA and IIB are related by T-duality. The link be- 
tween the field variables of the two models dictated by 
T-duality ^ enables one to derive a certain fractionally 
linear map ttba between their (Einstein-frame) Kasner 
exponents, which can be used, as above, to prove the 
Kasner-stability equivalence of the types IIA and IIB 
theories. Since type IIA is unstable, type IIB is also 
unstable. 

At this stage, we know that At, IIA and IIB are equiv- 
alent with respect to Kasner stability, and are all unsta- 
ble. It remains to tackle the type I and heterotic theo- 
ries, which are equivalent because their stability condi- 
tions are algebraically mapped onto each other by the 
5-duality transformation p 1 ^ = —p 1 ^, , p\ — p\. To 
study the Kasner-stability of the heterotic theory, we 
found very convenient to replace the Einstein-frame Kas- 
ner exponents (p^,Pi) by their string-frame counterparts 
(aj 1 ). The link between the two is (in d + 1 spacetime 
dimensions, see, e.g. j|) 



d-l-a 



(d - 1) ctj - a 
d-l-a 



(8) 



with a = (y]; cti) — 1 and i — 1, . . . , d. In terms of 
the a's the Kasner sphere S^ -1 is simply the usual unit 
sphere, J2i i a i) 2 = 1- I n our case, d = 9 and one should 
add a label "ft" to both the p's and the a's. In terms of 
the string-frame exponents it is found that the /i-stability 
conditions are equivalent to the simpler inequalities > 



and a? 



a£ < 1 (where i,j, k are all different) 



subjected to the constraints J2 t (a^) 2 = 1. It is easy 
to verify that these inequalities can never hold when the 
space dimension is d = 9. In that case, the closest one 
comes to satisfying the inequalities is the isotropic point 
cti = 1/3 for which the second inequality is saturated. 
This concludes our proof that the heterotic model (and 
therefore also the type I one) is Kasner unstable. Finally 
the two blocks of theories (At, A, B) and (/, h) are both 
Kasner unstable, though for different algebraic reasons. 



^Note that the fact that the IIA field variables depend on 
one less variable than the M-ones is unimportant. What is 
important is the map (Q) and the fact that we have taken into 
account in the stability criteria all possible dangerous terms 
in a generic solution. 
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Our results so far show that the generic solution of the 
low-energy string models can never reach a monotonic 
Kasner-like behaviour. Following the BKL approach 
one can go further and study the evolution near a cos- 
mological singularity as a sequence of Kasner-like "free 
flights" interrupted by "collisions" against the "potential 
walls" corresponding to the various stability- violating ex- 
ponents g, e or b. We have studied this problem Q and 
found the following universal "collision law" giving the 
Kasner exponents p 1 ^ of the Kasner epoch following a 
collision in terms of the old ones: 



9^ = 1 - 2 



(w ■ p) (w ■ u) 
{w ■ w) 



f- - 2 



(w ■ p) 
(w ■ w) 



(9) 



Here, p M stands for pP = p v and p 1 = Pi. The scalar 
products are computed with the metric occurring 
in the quadratic form (pa]), namely, Goo = 1, Got = 0, 



G, 



1, while the vector u (entering (|3b|)) has 



"covariant" components uo = 0, itj = 1. Finally, the 
"contravariant" vector characterizes the "wall" re- 
sponsible for the collision and is defined in such a way 
that the corresponding exponent (<?(p), e(p) or b(p)) reads 



w(p) = 



G 



[E.g., for the wall associated 



with the electric exponent e\Ti{p) = Pi + P2 + P3, Wfj, 
reads wq — 0, Wi = 1 for i = 1,2,3 and Wi = for 
i > 3.] The result (|) (which is a rescaled geometrical re- 
flection in the hyperplane w^p^ = 0) applies uniformly to 
all possible "walls" : gravitational, electric or magnetic, 
ft generalizes particular results derived by many authors 

Summarizing: In all string models, the generic solution 
near a cosmological singularity for the massless bosonic 
degrees of freedom exhibits BKL-type oscillations, i.e. 
a (formally infinite) alternation of Kasner-epochs. The 
primary sources of this BKL behaviour are (i) the pres- 
ence of p-forms in the field spectrum of the theories and, 
(ii) the strength of their dilaton couplings. In the ab- 
sence of p-forms, or if the A p 's were somewhat smaller, 
the monotonic Kasner behaviour would be stable and 
generic. The general rule defining the change of Kasner 
exponents from one epoch to the next is given by Eq. (^), 
where w is the "wall" (among the various gravitational, 
electric or magnetic ones) for which w(p) = w^p^ is most 
negative. We anticipate that the discrete dynamics (|]) 
will define (in all string models) a chaotic motion on the 
Kasner sphere. At this stage, the physical consequences 
of such a chaotic motion are unclear. It might consti- 
tute a problem for the pre-big-bang scenario Q which 
strongly relies on the existence, near a (future) cosmolog- 
ical singularity, of relatively large, quasi-uniform patches 
of space following a monotonic, dilaton-driven Kasner be- 
haviour. By contrast our findings suggest that the spatial 
inhomogeneity continuously increases toward a singular- 
ity, as all quasi-uniform patches of space get broken up 



into smaller and smaller ones by the chaotic oscillatory 
evolution. In other words, the spacetime structure tends 
to develop a kind of "turbulence" 113 ]. 

We are aware of the limitations of our result (tree- 
level bosonic massless modes only) but we think that 
our finding suggests that the full quantum, string-theory 
behaviour might be at least as complicated, near a cos- 
mological singularity, as our simplified analysis shows. 
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